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It is demonstrated that event generators LEPTO and PYTHIA can be modified to 
describe some azimuthal modulations. The comparisons of results obtained with 
modified LEPTO with existing data in the current fragmentation region of SIDIS 
are presented for Cahn and Sivers effects as well as the predictions for the target 
fragmentation region. The predictions for Sivers effect in the Drell— Yan process 
obtained with modified PYTHIA are also presented. The concept of hadronization 
function is discussed. 

1. Introduction 

The spin (in) dependent azimuthal asymmetries arising in high energy reac- 
tions allow us to study the dynamical effects related to spin and transverse 
momentum of partons in target nucleon and in hadronization. 

The Drell- Yan lepton pair production is the simplest process which does 
not include the hadronization dynamics and within QCD can be described 
using as a nonperturbative input only the parton distribution hmctions a : 

dohl+ha ^ + i-+x „ J2(f,/hJ q /H 3 + 1-2) ® da«+^ + +*-. (1) 
i 

The predictions for the Sivers 2 effect for this process recently have been pre- 
sented in 3,4 . Here it will be demonstrated that similar results are obtained 
using modified PYTHIA event generator. 

More nonperturbative inputs are needed to describe the semi-inclusive 
DIS (SIDIS) processes within the QCD formalism. Namely, for the particles 



a In the following the notations of 1 are used. 
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produced in the current fragmentation region (CFR) of SIDIS one needs to 
introduce fragmentation function, D q (z): 

da e P ^hx ^J2f g ® da lq ^ lq ® D h q . (2) 

For the particles produced in the target fragmentation region (TFR) an- 
other nonperturbative input — the fracture functions, M^ N (x,z), are 
needed: 

da e P ^ehx „ J2 da e ^ ® M>) N . (3) 
q 

In practice it is not easy to separate the two regions at moderate beam en- 
ergies and final hadronic state invariant masses . An alternative approach 
which is able to describe the particles production in the whole phase space is 
based on the LUND string fragmentation model and adopted in the Monte 
Carlo event generators 6 ' 7 . Here the SIDIS cross section can be represented 
as 8 

da^ lhx da 1 ^ ® H h q/N , (4) 

q 

where the hadronization function, H^ N , describs the particle production 
from the system formed by the struck quark and target remnant. 

In 9 ' 1 the role of parton intrinsic motion in SIDIS processes in CFR 
within QCD parton model has been considered at leading order; intrinsic 

is fully taken into account in quark distribution functions and in the 
elementary processes as well as the hadron transverse momentum, pj^ , with 
respect to fragmenting quark momentum. 

The average values of k± for quarks inside protons and p± for final 
hadrons inside the fragmenting quark jet where fixed by a comparison with 
data on Cahn effect 11 - the dependence of the unpolarized cross section 
on the azimuthal angle between the leptonic and the hadronic planes. The 
single spin asymmetry (SSA) A^*'~* s) recently observed by HERMES 12 
and COMPASS 13 Collaborations was successfully described by Sivers mech- 
anism. 

Here it will be demonstrated that both Cahn and Sivers effects can be 
implemented into Monte Carlo event generators. 

2. Including Cahn effect in LEPTO 

In the simplest case, corresponding to LO approximation of parton model, 
event generation in LEPTO proceeds in several steps: 
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(1) the hard scattering kinematics is generated, 

(2) the active quark inside the nucleon is chosen according to the quark 
density function f q (x,Q 2 ), 

(3) the transverse momentum of the final quark is simulated with Gaus- 
sian fcx and flat ip distributions. Note that the transverse momen- 
tum of the final final quark is equal to that of initial quark for 
leading order hard subprocess. 

(4) the string fragmentation machinery of JETSET program 7 is applied 
to form the final hadrons. 

The Cahn effect 11 is a kincmatical effect arising due to the presence 
of nonzero intrinsic transverse momentum of quarks in the nucleon. In 
the general case of non collinear kinematics Mandelstam variables depend 
on the quark transverse momentum and its azimuthal angle and at order 
0(k±/Q) one has 



Eq. (5) shows that the azimuthal angle of the final quark (and of the string's 
end associated with the struck quark) is now modulated with amplitude 
depending on y, Q and k±. 

This effect can be introduced in the LEPTD event generator at the step 
(3) of the event generation, when the transverse momentum and azimuthal 
angle of the scattered quark are generated. To do this the generation of the 
quark transverse momentum, k± , is left unchanged and then the azimuthal 
angle is generated according to Eq. (5). This leads to azimuthal modulation 
of the string axis. The momentum conservation means that the transverse 
momentum of the quark is balanced by that of the target remnant, which 
in turn means that the azimuthal angle of the target remnant ip qq = <p + 
it. Hence, one expects that the azimuthal angle of the hadrons in the 
target fragmentation region (TFR), xf < 0, will be modulated with a 
phase shifted by tt with respect to that in CFR. 

Data on azimuthal dependencies of SIDIS covering a large xf range 
have been obtained by the EMC Collaboration 14 for a beam energy of 
280 GeV. The xf dependence of (cos (j>h) / 'wi{y) , where wi(y) = (2 — 
y)s/l - yj (l + (1 - y) 2 ), obtained by using modified LEPTO for EMC kine- 
matics are presented in Fig. 1 together with data points from 14 (left panel). 
The simulations has been done with LO setting of LEPTO (LST(8)=0) and 
with values of the parameters describing intrinsic kr (PARL(3)=0.5) and 
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Figure 1. Left: The xp dependence of (cos iph)/ w l (y) f° r charged hadrons compared 
with EMC data. Right: Predictions of modified LEPTO for xp dependence of (cos <f>h) for 
different hadrons produced in 12 GeV unpolarized SIDIS process. 

fragmentation px (PARL(21)=0.45) as adopted in 1 . 

The predictions of modified LEPTD for (cos 4>h) of different hadron (tt + , 
ir~ , 7T° and p) produced in SfDIS on a proton target at future CEBAF 
12 GeV facility at JLab 15 are presented in Fig. 1 (right panel). One can 
see from Fig. 1 and that the predicted mean value of cos 4>h in the CFR is 
negative {cos<ph)cFR < 0, while in the TFR is positive {cos^^tfr > 0, 
as suggested by arguments based on transverse momentum conservation. 

3. Including Sivers effect in LEPTO 

The azimuthal modulation of the string transverse momentum in the pre- 
vious section was due to Cahn effect - the dependence of the non planar 
lepton-quark scattering cross section on the quark azimuth. The quark 
distribution, f q (x, fcj_) itself is independent of quark azimuthal angle. 

The situation is different when one considers SIDIS on a transversely 
polarized nucleon. Now a correlation between transverse momentum of 
quark, k± and target transverse polarization S T of type S T ■ [P x k±] is 
possible - the so called Sivers effect 2 . 

The unpolarized quark (and gluon) distributions inside a transversely 
polarized proton can be written as: 

/,/pt {x, k ± ) = f q/p {x, k ± ) + i A^/^/pT (x, k±) S T ■ (P x k±) . (6) 
Eq. (6) implies 

fq/pi {x, k±) + f q/p i (x, k ± ) = 2f q/p (x, k ± ) , 

/g/pT {X, fej.) - fq/pl {X, fej.) = A^/g/pT {x, S T ■ (P X fej_) , (7) 



February 2, 2008 4:23 Proceedings Trim Size: 9in x 6in 



kotzinian'comol 



5 

where f q / p {x,k±) is the unpolarized parton density and A N f q / p t (x, k±) 
is referred to as the Sivers function. Notice that, as requested by parity 
invariance, the scalar quantity S T ■ (P x fej_) singles out the polarization 
component perpendicular to the P — k± plane. For a proton moving along 
— z and a generic transverse polarization vector S T —\S T \ (cos (f>s, sin <j>s, 0) 
one has: 

S T ■ (P x fej.) = \S T \ sim>- <p s ) = \S T \ sm<j>siv , (8) 

where {ip — 4>s) — 4>siv is the Sivers angle. 

The Sivers function for each light quark flavor q — u,d are parameterized 
in the following factorized form 1 : 

A N f q/p ,(x,k 1 _) = 2Af q (x)h(k 1 _)f q/p (x,k ± ) , (9) 

where 

, (a + b "(K+M 

Af q (x) = N q x"«(l - xf« i^±^ , (10) 

a q o q 

Mfcx) = V2^^ e - fc - /M2 , (11) 

where N q , a q , b q and M (GeV/c) are parameters. Then Eq. (6)can be 
rewritten as 

f q / p i(x,kj_) = f q /p(x,k±)[l + \S T \J\f q (x)h(k±) sin<j) Slv }. (12) 

Again, the Sivers effect is incorporated into LEPTD at the stage 3) of 
the event generation in the same way as for the Cahn effect but now the 
azimuthal angle is generated according to Eq. (12). For simulations the 
following set of parameters compatible with those obtained in 1 ' 3 have been 
used: N u = N a = 0.5, N d = N 3 = -0.2, a q = 0.3, b q = 2 and M 2 = 0.36 
(GeV/c) 2 . 

In Fig. 2 the results of simulation for HERMES experimental conditions 
are compared with observed Sivers asymmetries 12 (left panel). 

Future facilities as Electron Ion Colliders or upgraded JLab will have 
larger kinematic coverage and will offer the possibility of studying the Sivers 
effect also with hadrons produced in the TFR. As an example, the simula- 
tions have been done for 12 GeV electron SIDIS of a proton target. The 
DIS cut Q 2 > 1 (GeV/c) 2 and W 2 > 4 GeV 2 and a cut on the produced 
hadron transverse momentum Pt > 0.05 GeV/c was imposed. The predic- 
tions for xf dependence for JLab kinematics is presented in Fig. 2 (right 
panel). The x and Pt dependencies in the TFR are presented Fig. 3. 
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Figure 3. Predicted dependence of ¥>s on x, left panel, (pt, right panel) for 

different hadrons produced in the TFR (xp < —0.1) of SIDIS of 12 GeV electrons off a 
transversely polarized proton target. 



In Fig. 4 the results for Sivers asymmetry in the Drell -YAn process 
obtained with modified PYTHIA generator are presented for two different 
energies: planned GSI pp collider with y/s = lAAGeV (left panel) and 
RHIC y/s = 200GeV (right panel). The results' 3 are similar to that obtained 
in 3 - 4 . 



b Note, that here the sign convention of 4 is adopted. 
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Figure 4. Predicted dependence of for Drell-Yan process on xp . Left: yfs = 14.4 
GeV, right: = 200 GeV 

4. Discussion and Conclusions 

The advantage of this MC based approach compared to standard QCD 
factorized approach is the full coverage of produced hadron phase space. 
Figs. 1 and 2 demonstrate that the modified LEPTO event generator well 
describs the data in the CFR both for Cahn and Sivers asymmetries. One 
can notice in Fig. 1 that the integrated experimental value of (cos 4>h) for 
charged hadrons in the CFR is not compensated by that in TFR. It seems 
improbable that this imbalance can be compensated by larger values of 
(cos 4>h) of neutral hadrons at xp ~ — 1. 

Note, that in present approach a possible modifications of hadronization 
in the case of polarized target have been ignored. In 8 it was shown that the 
hadronization functions in principle depends on polarization states of struck 
quark and target nucleon even for production of (pseudo) scalar or unpolar- 
ized particles. This dependence cannot be neglected at moderate energies 
and new nonperturbative input — the polarized hadronization functions, 
Air 1 ,^, are needed to describe the polarized SIDIS. The expression for the 
SIDIS helicity asymmetry is then looks as: 



Al(x,z,Q 2 ) = 



v-> 2 1 n2\vh I r>2Ul , A«(z,Q 2 )Aff£ /N (x,z,C 
Eg e g?0>£ )2 )H^ /N {x,Z,Q f 2 )[l + 



(13) 

Since the hadronization functions depend on target nucleon, active 
quark and produced hadron variables the new correlations as (a): S L ■ 
[p_L x fej_], (b): s L ■ [pj_ x k±], (c): [S T x p±] • [pj, x &?] etc are pos- 
sible. This correlations cannot be present separately in the distribution 
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and fragmentation functions. They will induce the azimuthal asymmetries 
in the SIDIS of (a): unpolarized lepton off longitudinally polarized target, 
(b): longitudinally polarized lepton off unpolarized target, (c): unpolarized 
lepton off transversely polarized target etc. 

The new high statistic measurements in both CFR and TFR of SIDIS 
will allow to check the predictions of the approach presented here and better 
understand the effects of the quark intrinsic transverse momentum and 
hadronization mechanism in SIDIS. The study of single spin asymmetries 
in Drell-Yan process will provide an additional test of our understanding 
od spin dependent phenomena. 
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